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Measurements of inclusive charged-hadron transverse-momentum and pseudorapid- 
^ ity distributions are presented for proton-proton collisions at = 0.9 and 2.36 TeV. 

The data were collected with the CMS detector during the LHC commissioning in 

December 2009. For non-single-diffractive interactions, the average charged-hadron 
^ transverse momentum is measured to be 0.46 ± 0.01 (stat.) ± 0.01 (syst.) GeV/c at 

. 0.9 TeV and 0.50 ± 0.01 (stat.) ± 0.01 (syst.) GeV/c at 2.36 TeV, for pseudorapidi- 

^ ties between —2.4 and +2.4. At these energies, the measured pseudorapidity den- 

Q sities in the central region, dNf,^/ df]\^^^^Q^, are 3.48 ± 0.02 (stat.) ± 0.13 (syst.) and 

^ 4.47 ± 0.04 (stat.) ± 0.16 (syst.), respectively. The results at 0.9 TeV are in agreement 

J> with previous measurements and confirm the expectation of near equal hadron pro- 

k>( duction in pp and pp collisions. The results at 2.36 TeV represent the highest-energy 

measurements at a particle collider to date. 



*See Appendix A for the list of collaboration members 
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1 Introduction 

Measurements of transverse-momentum (px) and pseudorapidity (//) distributions are reported 
for charged hadrons produced in proton-proton (pp) collisions at centre-of-mass energies {^/s) 
of 0.9 and 2.36 TeV at the CERN Large Hadron Collider (LHC) [IJ. The data were recorded 
with the Compact Muon Solenoid (CMS) experiment in December 2009 during two 2-hour 
periods of the LHC commissioning, demonstrating the readiness of CMS in the early phase of 
LHC operations. The results at -^s = 2.36 TeV represent the highest-energy measurements at a 
particle collider to date. 

The majority of pp collisions are soft, i.e., without any hard scattering of the partonic con- 
stituents of the proton. In contrast to the higher-px regime, well described by perturbative 
QCD, particle production in soft collisions is generally modelled phenomenologically to de- 
scribe the different pp scattering processes: elastic scattering, single-diffractive and double- 
diffractive dissociation, and inelastic non-diffractive scattering ||2l. 

The measurements presented in this paper are the inclusive primary charged-hadron multiplic- 
ity densities (dNjh/dpx arid dN^h/ drj) in the pseudorapidity range |?/| < 2.4, where px is the 
momentum of the particle transverse to the beam axis, and where N^^ is the number of charged 
hadrons in any given ?/ or px interval. The pseudorapidity f] is defined as — ln[tan(0/2)], where 
9 is the polar angle of the particle with respect to the anti-clockwise beam direction. 

Primary charged hadrons are defined as all charged hadrons produced in the interactions, in- 
cluding the products of strong and electromagnetic decays, but excluding products of weak 
decays and hadrons originating from secondary interactions. In this paper, the multiplicity 
densities are measured for inelastic non-single-diffractive (NSD) interactions to minimize the 
model dependence of the necessary corrections for the event selection, and to enable a com- 
parison with earlier experiments. The event selection was therefore designed to retain a large 
fraction of inelastic double-diffractive (DD) and non-diffractive (ND) events, while rejecting all 
elastic and most single-diffractive dissociation (SD) events. 

Measurements of dNch/ dpj and dN^h/ dr] distributions and their ^/s dependence are important 
for understanding the mechanisms of hadron production and the relative roles of soft and hard 
scattering contributions in the LHC energy regime. Furthermore, the measurements at the 
highest collision energy of 2.36 TeV are a first step towards understanding inclusive particle 
production at a new energy frontier. These measurements will be particularly relevant for the 
LHC as, when it is operated at design luminosity, rare signal events will be embedded in a 
background of more than 20 near-simultaneous minimum-bias collisions. These results will 
also serve as a reference in the measurement of nuclear-medium effects in PbPb collisions at 
the LHC. The differences in these distributions between pp and pp collisions are expected to 
be smaller than the attainable precision of these measurements ||3l- The results reported here at 
^/s = 0.9 TeV can therefore be directly compared to those previously obtained in pp collisions. 

This paper is organized as follows. In Section |2| the elements of the CMS detector relevant to 
this analysis are outlined. In Sections|3]and|4| the event selection and reconstruction algorithms 
are described. Results on dN^h/ dpj and dN^^^f drj are presented in Section|5]and compared with 
previous pp and pp measurements in Section |6] 

2 The CMS detector 

A detailed description of the CMS experiment can be found in Ref. [4J. The central feature 
of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a 
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uniform magnetic field of 3.8 T. Immersed in the magnetic field are the pixel tracker, the 
silicon-strip tracker (SST), the lead-tungstate crystal electromagnetic calorimeter (ECAL) and 
the brass /scintillator hadron calorimeter (HCAL). In addition to barrel and end-cap detectors 
for ECAL and HCAL, the steel/ quartz-fibre forward calorimeter (HE) covers the region of \f]\ 
between 2.9 and 5.2. The HE tower segmentation in t] and azimuthal angle cp (expressed in 
radians) is 0.175 x 0.175, except for \tj\ above 4.7 where the segmentation is 0.175x0.35. Muons 
are measured in gas-ionization detectors embedded in the steel return yoke. 

The tracker consists of 1440 silicon-pixel and 15 148 silicon-strip detector modules and mea- 
sures charged particle trajectories within the nominal pseudorapidity range \rj\ < 2.5. The 
pixel tracker consists of three 53.3 cm long barrel layers and two end-cap disks on each side of 
the barrel section. The innermost barrel layer has a radius of 4.4 cm, while for the second and 
third layers the radii are 7.3 cm and 10.2 cm, respectively. The tracker is designed to provide an 
impact-parameter resolution of about 100 fim and a transverse-momentum resolution of about 
0.7% for 1 GeV/c charged particles at normal incidence (//=0) O. 

During the data-taking period addressed by this analysis, 98.4% of the pixel and 97.2% of the 
SST channels were operational. The fraction of noisy pixel channels was less than 10^^. The 
signal-to-noise ratio in the SST depends on the sensor thickness and was measured to be be- 
tween 28 and 36, consistent with the design expectations and cosmic-ray measurements f5l,'6^. 
The tracker was aligned as described in Ref. [7J using cosmic ray data prior to the LHC com- 
missioning. The precision achieved for the positions of the detector modules with respect to 
particle trajectories is 3-4 }im in the barrel for the coordinate in the bending plane. 

Two elements of the CMS detector monitoring system, the Beam Scintillator Counters (BSCs) lH 
m and the Beam Pick-up Timing for the experiments (BPTX) devices [4, 9J, were used to trigger 
the detector readout. The two BSCs are located at a distance of ±10.86 m from the nominal 
interaction point (IP) and are sensitive in the \t]\ range from 3.23 to 4.65. Each BSC is a set of 
16 scintillator tiles. The BSC elements have a time resolution of 3 ns and an average minimum- 
ionizing-particle detection efficiency of 96.3%, and are designed to provide hit and coincidence 
rates. The two BPTX devices, located around the beam pipe at a distance of ±175 m from the IP 
on either side, are designed to provide precise information on the bunch structure and timing 
of the incoming beam, with better than 0.2 ns time resolution. 

The CMS experiment uses a right-handed coordinate system, with the origin at the nominal 
interaction point, the x axis pointing to the centre of the LHC, the y axis pointing up (perpen- 
dicular to the LHC plane) and the z axis along the anticlockwise-beam direction. The azimuthal 
angle, (p, is measured in the {x,y) plane, where ^ = is the +x and (p = n/2\s the +y direction. 

The detailed Monte Carlo (MC) simulation of the CMS detector response is based on Geant4 
IITOl . The position and width of the beam spot in the simulation were adjusted to that deter- 
mined from the data. Simulated events were processed and reconstructed in the same manner 
as collision data. 

3 Event selection 

This analysis uses two LHC collision data sets collected with pp interaction rates of about 11 
and 3 Hz at a/s = 0.9 and 2.36 TeV, respectively. At these rates, the probability for more than 
one inelastic collision to occur in the same proton bunch crossing is less than 2 x 10^^ at both 
collision energies. 

Events were selected by a trigger signal in any of the BSC scintillators, coincident with a sig- 
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nal from either of the two BPTX detectors indicating the presence of at least one proton bunch 
crossing the IP. From these samples, collision events were selected offline by requiring BPTX 
signals from both beams passing the IP and at least one reconstructed charged particle trajec- 
tory in the pixel detector originating from within 0.2 cm of the beam position in the transverse 
direction (Section [4. The total number of collision events and the numbers of collision events 
passing each requirement are listed in Table [T] 

Table 1: Numbers of events per data sample used in this analysis. The offline event selection 
criteria are applied in sequence, i.e., each line includes the selection of the lines above. 



Centre-of-mass Energy 


0.9 TeV 


2.36 TeV 


Selection 


Number of Events 


BPTX Coincidence + one BSC Signal 


72637 


18 074 


One Pixel Track 


51308 


13 029 


HF Coincidence 


40 781 


10 948 


Beam Halo Rejection 


40 741 


10 939 


Beam Background Rejection 


40647 


10 905 


Valid Event Vertex 


40320 


10 837 



To select NSD events, a coincidence of at least one HF calorimeter tower with more than 3 GeV 
total energy on each of the positive and negative sides of the HF was required. Events con- 
taining beam-halo muons crossing the detector were identified by requiring the time difference 
between any two hits from the BSC stations on opposite sides of the IP to be within 73 ± 20 ns. 
Such events were removed from the data sample. Beam-induced background events produc- 
ing an anomalously large number of pixel hits were excluded by rejecting events with pixel 



clusters (Section 4.2 1 inconsistent with a pp collision vertex. This rejection algorithm was only 
applied for events with more than 150 pixel clusters, providing a clean separation between col- 
lision events and beam background events. Finally, events were required to contain at least one 
reconstructed primary vertex, as described in Section |4] 

To study beam-induced background, the event selection criteria were also applied to a data 
sample obtained by selecting events with only a single unpaired bunch crossing the IP. The 
contamination of background events in the colliding-bunch data sample was estimated by tak- 
ing into account the total unpaired and paired bunch intensities and was found to be negligible 
(<0.1%). The total number of cosmic-ray muons in the selected data sample was estimated to 
be less than one event, and was also neglected. 

The event selection criteria are expected to have high efficiency for the NSD part of the pp cross 
section, while rejecting a large fraction of the SD component of pp interactions. The efficiency 
of the event selection for the different processes and centre-of-mass energies was determined 
using simulated events obtained from the pythia IITTl (version 6.420, tune D6T, l|T2l ) and PHO- 
JET HTSl [T4| (version 1.12-35) event generators processed with a MC simulation of the CMS 
detector response (hereafter simply called pythia and phojet). In the case of phojet, the dis- 
cussion and numerical values concerning the DD process given in this paper contain both the 
DD and the Double-Pomeron-Exchange (DPE) processes. The relative event fractions of SD, 
DD and ND processes and event selection efficiencies at = 0.9 and 2.36 TeV are listed in 
Table |2]for these two samples. 

The measurements were corrected for the selection efficiency of NSD processes and for the 
fraction of SD events contained in the data sample after the event selection. Based on the 
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PYTHIA (phojet) event generator, the fractions of SD events contained in the selected data 
samples were estimated to be 5.2% (4.9%) at 0.9 TeV and 6.3% (5.0%) at 2.36 TeV. 

Table 2: Expected fractions of SD, DD, ND and NSD processes ("Frac") obtained from the 
PYTHIA and PHOJET event generators before any selection and the corresponding selection ef- 
ficiencies ("Sel. Eff.") determined from the MC simulation. 



PYTHIA 


PHOJET 


Energy 


0.9 TeV 


2.36 TeV 


0.9 TeV 


2.36 TeV 


Frac. 


Sel. Eff. 


Frac. 


Sel. Eff. 


Frac. 


Sel. Eff. 


Frac. 


Sel. Eff. 


SD 


22.5% 


16.1% 


21.0% 


21.8% 


18.9% 


20.1% 


16.2% 


25.1% 


DD 


12.3% 


35.0% 


12.8% 


33.8% 


8.4% 


53.8% 


7.3% 


50.0% 


ND 


65.2% 


95.2% 


66.2% 


96.4% 


72.7% 


94.7% 


76.5% 


96.5% 


NSD 


77.5% 


85.6% 


79.0% 


86.2% 


81.1% 


90.5% 


83.8% 


92.4% 



The generated charged-hadron multiplicity distribution is shown in Fig.jlp in the range \ < 2.5 
for all inelastic events after event selection. The event selection efficiency for NSD events is 
shown in Fig.jlJ? as a function of generated charged-hadron multiplicity in the region \f]\ < 2.5. 
The correction for the event selection efficiency was applied as a function of number of recon- 
structed charged particles per event, as illustrated at generator level in Figs.jl^ andjl]?. 

The sum of the corrections to the dN^h/ drj measurements due to the NSD event selection ef- 
ficiency and the SD event contamination typically amounts to 8%. The corrections applied in 
the analysis are based on PYTHIA using the default SD and DD process fractions as listed in 
Tabled 

The PYTHIA predictions for the SD and DD fractions differ from those of PHOJET, and are not 
fully consistent with existing measurements, as explained in Section|5j These differences propa- 
gate to a systematic uncertainty of 2% in the dNch/ drj measurement. To estimate the additional 
systematic uncertainty on the event selection efficiency correction resulting from the possible 
inaccuracies in the detector simulation, the analysis was repeated after replacing the HE event- 
selection criterion with a two-sided hit coincidence of signals in the BSC detectors. Based on 
this comparison, an additional 1% systematic uncertainty was assigned to the dNf^^/drj mea- 
surements. 



4 Reconstruction algorithms 

The analysis presented in this paper measures the dN^h/ dtj and dJVch/ dpj distributions of pri- 
mary charged hadrons. The dN^^^f drj distributions were obtained using three methods based 
on counting of (i) reconstructed clusters in the pixel barrel detector; (ii) pixel tracklets com- 
posed of pairs of clusters in different pixel barrel layers; and (iii) tracks reconstructed in the full 
tracker volume, combining the pixel and strip hits. The cluster counting method provides an 
independent measurement for each pixel barrel layer, and the tracklet method for each pair of 
layers. The third method also allows a measurement of the dN^-h/dpj distribution. All three 



methods rely on the reconstruction of the primary vertex (PV) described in Section 4.1 



The pixel-cluster-counting method has the advantage of having the largest pj acceptance down 
to small transverse momentum (30 MeV/c), is insensitive to geometrical misalignment of the 
detector and does not require detailed knowledge of the primary vertex position resolution. 
A potential disadvantage is the sensitivity to backgrounds from collisions with residual gas in 



4.1 Primary vertex reconstruction 
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Figure 1: (a) Generated multiplicity distributions of primary charged hadrons in the range 
\f]\ < 2.5 for ^ys =0.9 TeV (solid dots and histogram) and 2.36 TeV (open circles and dashed 
histogram) after the event selection is applied to the reconstructed events, using inelastic events 
from the PYTHIA (histograms) and PHOJET (symbols) event generators, (b) The event selection 
efficiency expected for NSD events from the PYTHIA (histograms) and PHOJET (symbols) event 
generators as a function of generated charged hadron multiplicity in the region \fl\ < 2.5. 

the beam pipe (beam-gas collisions), from secondary particles produced in the detector mate- 
rial and from low-px particles curling in the axial magnetic field (loopers). The pixel- tracklet 
method is capable of measuring and correcting for the combinatorial background and has a pj 
threshold of 50 MeV/c. The third method uses the tracker (pixel and SST) to build tracks. It 
requires at least two pixel hits in different layers, has the largest pj threshold (~100 MeV/c) 
and algorithmic complexity, but is the most robust against background hits produced by par- 
ticles not originating from the collision. The charged-particle multiplicity was corrected in all 
three methods for the small contamination (< 1%) of primary charged leptons. The measured 
dNch/ dr] values were evaluated by extrapolating or correcting to px = for all the three analy- 
sis methods. 

The three reconstruction methods are described in Sections 14.21 - 14.41 
4.1 Primary vertex reconstruction 

The X, y and z positions of the luminous region where protons of both beams interact, hereafter 
referred to as beam spot, are obtained for each data set from three-dimensional vertex fits based 
on tracks reconstructed with pj > 0.9 GeV/ c, using the full event sample. The RMS of the beam 
spot in the transverse directions was found to be less than 0.05 cm. The beam spot position and 
dimensions were found to be stable within a given data set. 

To reconstruct the z coordinate of the PV for each event, tracks consisting of triplets of pixel 
hits were formed. The minimum transverse momentum of these tracks is « 75 MeV/c. The 
tracks were required to originate from the vicinity of the beam spot with a transverse impact 
parameter (dj) smaller than 0.2 cm. Of these tracks, only those with dj < A aj, where cfj is 
the quadratic sum of the uncertainty in dj and the RMS of the beam spot in the transverse 
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Zpv [cm] Angle-corrected cluster charge [ke] 



Figure 2: |(a)| The distribution of the reconstructed z position of the primary vertex in the data 
(dots), compared to that from the PYTHIA simulation (histogram). |(b)| Distribution of the cluster 
charge multiplied by | sin 9 \ in the data (dots) and the simulation (histogram), for the clusters 
selected for analysis. 

direction, were used in the vertex reconstruction. 

The vertex-reconstruction algorithm uses the z coordinate of the tracks at the point of closest ap- 
proach to the beam axis and the corresponding estimated measurement uncertainty (cTz). It per- 
forms an agglomerative clustering by adding tracks to form groups. These groups (denoted the 
and;'* group) are then merged based on their normalized distance, d?- = (z; — zy)^/ {af + aj) 
where a, and dj are the uncertainties of the z,- and Zy positions, with a fast nearest-neighbour 
search algorithm [15]. The z position and its uncertainty (Tz for the newly joined group are cal- 
culated using a weighted average. The clustering process stops when the smallest normalized 
distance between the remaining groups gets larger than 8.0, where the stopping condition was 
optimized using simulated events. Only vertices formed from at least two tracks were con- 
sidered except when only one track was reconstructed in the event. In the latter case the PV 
position was defined as the point of the closest approach of the track to the beam axis. The 
fraction of single-track vertices in the selected data sample is 1.7% (1.3%) at 0.9 TeV (2.36 TeV). 
The overall PV reconstruction efficiency, evaluated from the data after all other event selection 
cuts, is m excess of 99% and the fraction of events with more than one primary vertex candi- 
date is 5.0% (7.4%) at 0.9 TeV (2.36 TeV). In the rare case of multiple PV candidates, the vertex 
composed of the largest set of tracks was chosen. 

The reconstructed primary vertex resolution in the z direction is a function of the associated 
track multiplicity (N) and was found to be parameterized adequately as 0.087 cm / N^-^ using 
simulated events. 

The distribution of the reconstructed z position of the PV is shown in Fig.|2^. Overlaid is the PV 
distribution in simulated events, the position and RMS of which were adjusted to reproduce 
the beam spot measured in data. 



4.2 Pixel cluster counting method 



7 



4.2 Pixel cluster counting method 

The pseudorapidity distribution of primary charged hadroris produced in a pp collision can 
be measured by counting the number of clusters they create when traversing each of the three 
pixel barrel layers and applying appropriate corrections, as described in this section. 

The energy deposited by charged particles traversing a pixel detector layer is spread over mul- 
tiple pixels. Adjacent pixels with a charge measurement above a readout threshold of typically 
2740 electrons are combined into pixel clusters to integrate the total charge deposit [6J. The 
cluster size and charge depend on the incident angle of the particle with respect to the active 
detector surface. The cluster length along the z axis ranges from 1-2 pixels at normal incident 
angle up to 14-16 pixels at shallow crossing angles. Figure |2]3 shows the measured distribution 
of cluster charge multiplied by | sin 9\ (or 1 / cosh rj) after the cluster selection discussed below, 
compared to the simulation. Here, 6 is the polar angle of the straight line connecting the PV to 
the cluster. 

The peak position is consistent with the expected charge of 22 ke, while the width of the 
distribution is slightly larger in the data than in the simulation due to gain calibration non- 
uniformities. 

The cluster counting method correlates the observed pixel-cluster length in the z direction, 
expressed in number of pixels, with the expected path length traveled by a primary particle at a 
given rj value. For primary particles the cluster length in z is proportional to | cot 6\ (or | sinh J/|) 
as displayed in Fig. |3^. Small clusters at large \ are due to loopers, secondary particles and 
daughters of long-lived hadrons. Clusters from these background particles were efficiently 
removed by the cluster-length cut represented by the solid line in Fig. |3^. To allow for an 
efficient background rejection, only the barrel part of the pixel detector was used, where the 
detector units are parallel to the beam axis, as opposed to the pixel end-caps. Furthermore, the 
rj range for the cluster counting was restricted to < 2 to avoid acceptance problems due to 
the slightly off-centred position of the luminous region. 

The event selection efficiency and the SD contribution for a given total multiplicity of selected 
clusters (M) for each pixel barrel layer can be determined from Monte Carlo simulation. The 
overall change of the dN^.^/ drj value due to this correction is 9% for both collision energies. 

The fraction of clusters created by loopers above the cluster-length cut (1% and 5% for t] = 2 
and rj = 0.5, respectively) can be estimated by measuring the total number of clusters below 
the cut in data, corrected by the ratio of looper clusters below and above the cut in simulated 
events. The number of clusters eliminated by the cluster-length cut was found to be higher 
in data than in simulated events by 10-20% due to a slightly larger abundance of secondary 
particles and loopers, while the observed number and length distributions of clusters above 
the cut was found to agree with the simulation. 

The corrections for absorption in the beam pipe and detector material, secondary particles, 
daughters of long-lived hadrons, delta-ray electrons and double hits caused by geometrically 
overlapping detector units were all evaluated, in bins of and M, with simulated data. The size 
of these corrections is 10%, 23% and 41% for the first, second and third detector layer, respec- 
tively. Varying the charged-particle multiplicity in the event generator by 50% only causes a 
±3% relative change in these corrections. Their dependence on r] and pixel-cluster multiplicity 
is similarly small. 
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At] 



Figure 3: (a) Pixel cluster length along z as a function of t] for the 900 GeV data. The solid line 



illustrates the cut applied in the cluster counting method, (b) The Arj distribution of clusters on 
tracklets in the data (dots and circles), together with the distribution obtained from the pythia 
simulation (solid and dotted lines), for both 0.9 and 2.36 TeV collision energy. The dashed line 
shows the A?/ distribution of clusters for primary charged-particle tracks in the Monte Carlo 
simulation at 0.9 TeV. The tail of the At] distribution comes from the combinatorial background. 

4.3 Pixel-tracklet method 

This method was first used to measure charged-hadron multiplicities by the PHOBOS experi- 
ment at RHIC [16J. Pixel tracklets are constructed from combinations of two pixel hits in any 
two pixel barrel layers. The difference in the angular positions of the two clusters with respect 
to the PV, Af/ and A^, are calculated for each tracklet. If two tracklets share a hit, the track- 
let with the larger Arj is discarded. The At] distribution of reconstructed tracklets is shown in 
Fig. |3j3, together with the corresponding distribution from simulated data and a separate dis- 
tribution for simulated primary particles only. Tracklets from primary particles display a sharp 
peak at At] = 0, while the tracklets from the combinatorial background have an extended tail. 
The simulated At] distributions are in good agreement with data. 

To suppress the combinatorial background, only tracklets with \Arj\ < 0.1 and \A(p\ < 1.0 
were selected. Since the combinatorial background is flat in A^, the remaining fraction of 
background tracklets in the signal region \A(p\ < 1.0 can be estimated from tracklets with 
1.0 < I A^l < 2.0. This data-driven estimate of the background accurately describes the raw At] 
distribution of tracklets for \At]\ > 2, where no signal from primary particles is expected from 
the MC simulation. Typical values of this estimated background fraction in the signal region 
increase with \r] \ from 2% to 30%. The tj range for the tracklet method was restricted to \ f]\ < 2 
to avoid a large acceptance correction. 

The contribution from secondary particles, reconstruction efficiency and geometrical accep- 
tance needs to be accounted for to determine the number of primary charged hadrons. These 
correction factors were calculated using PYTHIA simulations for background-subtracted track- 
lets in bins of z position of the PV, pseudorapidity, and tracklet multiplicity. The magnitude of 
the correction varies with 1 1] \ from to 20%. The correction factors were also cross-checked by 
PHOJET simulations and only cause a 2-3% change in the dN^h/ dtj result. 



4.4 Tracking method 
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The correction for the event selection efficiency and the SD contribution was determined for 
each tracklet multiplicity bin. The overall change in the dN^h/ drj value due to this correction is 
about 8% at ?/ = 0. 

4.4 Tracking method 

Pixel and SST detectors were used to reconstruct tracks, including both barrel and end-cap lay- 
ers. The acceptance was limited to |//| < 2.4 to avoid edge effects. The iterative reconstruction 
procedure described below follows Refs. (iTl [TSl , but was further optimized for primary-track 
reconstruction in minimum bias events. 

In the first step of track reconstruction, tracks with three pixel hits (triplets) are built using the 
X and y positions of the beam spot and the z coordinate of the primary vertex as constraints. 
These clean pixel tracks are used as seeds for the Kalman-filter-based trajectory-building al- 
gorithm in the SST. The resulting trajectories are stored. Before the second tracking step, the 
pixel and strip hits associated with the tracks found in the first step are removed from further 
consideration. The second step uses pixel triplet seeds as well, but does not require a vertex 
constraint and has a looser transverse impact parameter requirement than in the first step. Af- 
ter removal of hits associated with tracks found in the second step, the third tracking step finds 
primary tracks seeded by two hits in the pixel detector. At least three hits were required for a 
track to be accepted. 

Tracks found during the three iterative steps were collected and a second iteration of the PV 



reconstruction, as described in Section 4.1 was performed to refine primary vertex position 
determination. Finally, the tracks were refit with the corresponding vertex constraint, thus 
improving their t] and pj resolution. 

In this analysis, a reconstructed track was considered as a primary-track candidate if it is com- 
patible with originating from the PV (dj < min(4t7x, 0.2 cm) and dz < 4(7z, where dz is the 
distance between the point of the closest approach of the track to the beam axis and the PV 
along the z direction). 

Studies with simulated events showed that the combined geometrical acceptance and recon- 
struction efficiency for the tracking method exceeds 50% around pj ~ 0.1, 0.2 and 0.3 GeV/c 
for pions, kaons and protons, respectively. The efficiency is about 96% in the | ;y | < 1 region for 
Pj > 0.25 GeV/c, and is above 80% for pions at pj = 0.15 GeV/c. By requiring the geometrical 
shapes of the pixel clusters to be consistent with the crossing angle and direction of the track, 
the fraction of fake tracks was kept below 1%. The fraction of duplicated tracks (e.g., from 
loopers) was estimated to be about 0.1% in the central region, rising to 0.5% at large |7/|. 

The measured yield in data was corrected, based on MC simulation and comparisons with 
data, for geometrical acceptance (2% correction for pj > 200 MeV/ c), efficiency of the recon- 
struction algorithm (5-10% for pj > 300 MeV/c), fake and duplicate tracks (<1% each). The 
contamination of less than 2% from decay products of long-lived hadrons, photon conversions 
and inelastic hadronic interactions with the detector material was also subtracted. To obtain 
the dNch/ df] result from the pj spectrum, an extrapolation to px = was necessary, resulting 
in an increase of 5% in the estimated number of charged hadrons. 

Corrections based on the average hit efficiency of pixel layers, size of the beam spot, longitudi- 
nal and transverse impact-parameter resolutions of pixel tracks were validated with data. As 
an example, the average number of pixel and strip hits found on tracks in the range \ri\ < 1 is 
shown in Figs. |4^ and |^, together with the expectation from PYTHIA. Somewhat fewer par- 
ticles are predicted with pj < 500 MeV/ c than seen in the data, which results in the small 
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Figure 4: 



(a) The distribution of the number of pixel hits attached to reconstructed tracks in 
< 1 (closed circles), compared to the CMS detector simulation (histogram), (b) 



the region of | rj 

The distribution of the number of hits in the SST detector attached to reconstructed tracks in 
the region oi\r]\ < 1 (closed circles), compared to the CMS detector simulation (histogram). 

difference in the number of tracks with few SST hits in Fig. |^. This small difference, which 
originates from limitations of the PYTHIA generator, does not affect the final measurement. The 
correction for the event selection efficiency and the SD contribution was determined for each 
track multiplicity bin, and has an overall magnitude of 8.3%. 



5 Results 

5.1 Systematic uncertainties 

Various corrections and their event-selection and model dependence contribute to the system- 
atic uncertainties of the measured quantities. A summary of these systematic uncertainties 
averaged over rj and pj is given in Table |3] and discussed below. 

The uncertainties related to the trigger bias and to the event selection are common to all the 
analysis methods. The efficiency of the trigger and event-selection was corrected for by the 
prediction of the PYTHIA D6T event generator combined with the full Geant4 simulation of 
the CMS detector. The material description relevant for this analysis was verified by studies 
of photon conversion probabilities in the data, found to be in agreement with those obtained 
from the simulation. 

Because single- and double-diffractive pp collisions have much smaller charged-hadron mul- 
tiplicities per event than non-diffractive events, they contribute to the uncertainty in the mea- 
sured dNch/dt] mostly through the uncertainty in the fraction of SD and DD events passing 
the event selection criteria. The fractions of SD events for ^/s = 0.9 TeV in PYTHIA and PHO- 
JET (Table |2]) are 23% and 19%, respectively. The UA5 experiment measured 15.5% for this 
fraction 11191 . Based on the simulated trigger efficiencies for the different event types, only 
5.5% of events passing the analysis event selection are expected to be single-diffractive events. 
From the aforementioned variations of SD fractions, an uncertainty of ±1% is attributed to this 
correction. The contribution of the uncertainty of the fraction of DD events was estimated sim- 



5.2 Charged hadron transverse-momentum distributions 
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ilarly to be ±1%. Since underestimated DD and SD fractions both lead to an underestimated 
dNch/ dt] result, a conservative linear sum of 2% was assigned to the above systematic uncer- 
tainty. The trigger efficiency of the BSC is more than 98% for events with a valid vertex, and 
even a 5% uncertainty in single-particle detection efficiency of its individual segments results 
in a negligible uncertainty in the final result. The trigger efficiency of the BSC and the event se- 
lection efficiency of the HF detector were both measured from data and found to be consistent 
within 1% with the MC simulation. The total systematic uncertainty from propagating all event 
selection and trigger related uncertainties is 3%. The measurement of the average transverse 
momentum is less sensitive to the trigger selection efficiency. A smaller, 1% uncertainty was 
therefore assigned to that result. 

The geometrical acceptance was studied by comparing the hit occupancy of the pixel barrel 
with the predictions from the simulation. The efficiency of the pixel hit reconstruction was 
estimated using tracks propagated from the SST to the pixel detector and by extrapolating pixel 
tracklets to the unused pixel barrel layer. The measured pixel hit efficiency was found to exceed 
99% with a 0.5% uncertainty from both methods, which propagates into 0.5% uncertainty in 
the pixel-counting-based, 1% in the tracklet-based, and 0.3% in the track-based results. If the 
collected charge in one or more pixels in a cluster remains below the threshold, the cluster may 
be split. The splitting rate was estimated from the geometrical distance distributions of close- 
by pixel clusters found in the data and in the Monte Carlo simulation and found to be 0.5-0.9% 
in the simulation and 1.0-1.5% in data. 

The uncertainty related to the cluster and tracklet selections was estimated by varying the se- 
lection cuts. An additional 3% and 2% uncertainty was assigned to the tracklet and track recon- 
struction algorithm efficiencies, respectively. Corrections for loopers and secondary particles 
are simulation dependent; the tracklet- and pixel-counting-based methods have low rejection 
power compared to the tracking method, thus carry a larger systematic uncertainty (as shown 



The effects of the geometrical misalignment of the pixel barrel detector were simulated and a 
1% uncertainty was assigned to the results from the tracklet-based method. Hits from beam- 
induced backgrounds coinciding with the collision were estimated to be very rare, and a con- 
servative 1% random hit contribution was propagated to obtain the uncertainty of the results. 
The corrections for multiple track counting and fake track rate were estimated from the Monte 
Carlo simulation and found to be less than 1%. The uncertainty of the extrapolation to the full 
pT-range depends on the low-pj reach of the three methods and varies between 0.2 and 0.5%. 
While the sources of uncertainties are largely independent from each other, they are correlated 
among all the data points. 

5.2 Charged hadron transverse-momentum distributions 

Tracks with \ f]\ < 2.4 and pj > 0.1 GeV/c were used for the measurement of dN(,\y/dpj. The 
measured average charged-hadron yields per NSD event are shown in Fig. |5^, as a function of 
Pt in bins of \rj\. The yields were fit by the Tsallis function (Eq. [l]|, which empirically describes 
both the low-px exponential and the high-px power-law behaviours [20.. 2H : 



in Table |3]|. 



dp^ 



'ch 



iTifT p dtjdpj 



1 £ d'^N^^ 




(1) 



where y = 0.5 ln[(£ + p^)/ {E — pz)] is the rapidity; C(n, T, m) is a normalization constant that 
depends on n, T and m; Ej = + pj — m and m is the charged pion mass. This function 
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5 Results 



Table 3: Summary of systematic uncertainties. While the various sources of uncertainties are 
largely independent, most of the uncertainties are correlated between data points and between 
the analysis methods. The event selection and acceptance uncertainty is common to the three 
methods and affects them in the same way. The values in parentheses apply to the {pr) mea- 
surement. 



Source 


Pixel Counting [%] 


Tracklet [%] 


Tracking [%] 


Correction on event selection 


3.0 


3.0 


3.0 (1.0) 


Acceptance uncertainty 


1.0 


1.0 


1.0 


Pixel hit efficiency 


0.5 


1.0 


0.3 


Pixel cluster splitting 


1.0 


0.4 


0.2 


Tracklet and cluster selection 


3.0 


0.5 




Efficiency of the reconstruction 




3.0 


2.0 


Correction of looper hits 


2.0 


1.0 




Correction of secondary particles 


2.0 


1.0 


1.0 


Misalignment, different scenarios 




1.0 


0.1 


Random hits from beam halo 


1.0 


0.2 


0.1 


Multiple track counting 






0.1 


Fake track rate 






0.5 


Pt extrapolation 


0.2 


0.3 


0.5 


Total, excl. common uncertainties 
Total, incl. common uncert. of 3.2% 


4.4 
5.4 


3.7 
4.9 


2.4 
4.0 (2.8) 



provides both the inverse slope parameter T, characteristic for low pj, and the exponent n, 
which parameterizes the high-px power-law tail. These fit parameters change by less than 5% 
with f], thus a fit to the whole region \t]\ < 2.4 was performed. The px spectrum of charged 
hadrons, 1/ {2npj)d^Nch/ dtjdpj, in the region \ f]\ < 2.4, was also fit with the empirical function 
(Eq. [l]| and is shown in Fig. |^. The px resolution of the CMS tracker was found to have a 
negligible effect on the measured spectral shape and was therefore ignored in the fit function. 
For the 0.9 TeV data, the inverse slope parameter and the exponent were found to be T = 0.13 ± 
0.01 GeV and n = 7.7 ± 0.2. For the 2.36 TeV data, the values were T = 0.14 ± 0.01 GeV and n = 
6.7 ± 0.2. The average transverse momentum, calculated from the measured data points adding 
the low- and high-px extrapolations from the fit is (px) = 0.46 ± 0.01 (stat.) ± 0.01 (syst.) GeV/c 
for the 0.9 TeV and 0.50 ± 0.01 (stat.) ± 0.01 (syst.) GeV/c for the 2.36 TeV data. 

The dNch/df] spectrum was obtained by summing the measured differential yields for 0.1 < 
Px < 3.5 GeV/c and adding the result to the integral of the fit function for px < 0.1 GeV/c and 
Px > 3.5 GeV/c. The latter term amounts to 5% of the total. 

5.3 Charged hadron pseudorapidity density 

The summary of results on the pseudorapidity density distribution of charged hadrons is 
shown in Fig.|6] The dNch/ dtj results for the three layers in the cluster-counting method and the 
three layer-pairs in the pixel-tracklet method are consistent within 3%. These results from the 
various layers and from the different layer pairs were combined to provide one set of data from 
each analysis method, as shown in Fig.|6^. The error bars include the systematic uncertainties 
of about 2.4-4.4% specific to each method, estimated from the variations of model parame- 
ters in the simulation used for corrections and the uncertainties in the data-driven corrections. 
The systematic uncertainties common to all the three methods, which amount to 3.2%, are not 
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Figure 5: (a) Measured differential yield of charged hadrons in the range \r]\ < 2.4 in 0.2- 
unit-wide bins of \t]\ for the 2.36 TeV data. The measured values with systematic uncertainties 
(symbols) and the fit functions (Eq. are displayed. The values with increasing ?] are succes- 
sively shifted by four units along the vertical axis, (b) Measured yield of charged hadrons for 
|?/| < 2.4 with systematic uncertainties (symbols), fit with the empirical function (Eq.jl]). 

shown. The results from the three analysis methods are in agreement. The larger fraction of 
background hits in the data compared to simulation affects the cluster-counting method differ- 
ently from the other two, which results in a small difference at high ?], well accounted for by 
the systematic uncertainty of the measurement. 



6 Discussion 

The average transverse-momentum and pseudorapidity densities of charged hadrons derived 
from the measured data can be compared to results from earlier experiments as a function 
of the collision energy. The average transverse momentum of charged hadrons was obtained 
from the fits (Eq.jl]) to the transverse-momentum spectrum (Fig. [Sj?). At low energies the en- 
ergy dependence of (pj) can be described by a quadratic function of Ins. The (pj) from this 
measurement, shown in Fig. [7^, follows the general trend. At 0.9 TeV it is similar to the results 
from pp collisions at the same energy | |22| . 

The dNi,h/dtj distribution was calculated as the weighted average of the data from the three 
reconstruction methods, taking into account their systematic uncertainties, excluding the com- 
mon ones, as listed in Table |3] The averaged result is shown in Fig. ^ and compared to mea- 
surements at the same accelerator (ALICE, pp [23]) and to previous measurements at the same 
energy but with different colliding particles (UA5, pp [,24J ). The shaded error band on the CMS 
data indicates systematic uncertainties, while the error bars on the data from UA5 and ALICE 
display statistical uncertainties only. No significant difference is observed between the dN^h/ dtj 
distributions measured in pp and pp collisions at i/s = 0.9 TeV. 

The dNch/ dtj distribution is weakly //-dependent, with a slow increase towards higher rj val- 
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Figure 6: (a) Reconstructed dNch/ drj distributions obtained from the cluster counting (dots with 
error bars), tracklet (squares) and tracking (triangles) methods, in pp collisions at 0.9 TeV (filled 
symbols) and 2.36 TeV (open symbols). The error bars include systematic uncertainties (as dis- 
cussed in the text), excluding those common to all the methods, (b) Reconstructed dN^h/ drj 
distributions averaged over the cluster counting, tracklet and tracking methods (circles), com- 
pared to data from the UA5 [24| (open squares) and from the ALICE Il23l (open triangles) ex- 
periments at 0.9 TeV, and the averaged result over the three methods at 2.36 TeV (open circles). 
The CMS and UA5 data points are symmetrized in //. The shaded band represents systematic 
uncertainties of this measurement, which are largely correlated point-to-point. The error bars 
on the UA5 and ALICE data points are statistical only. 



ues, and an indication of a decrease at > 2 for the 0.9 TeV data. At 2.36 TeV, the entire 
distribution is wider due to the increased collision energy hence the larger t] range available 
for inclusive particle production. For \rj\ < 0.5, the corrected results average to dNch/ d)] = 
3.48 ± 0.02 (stat.) ± 0.13 (syst.) and dNa,/dt] = AA7 ± 0.04 (stat.) ± 0.16 (syst.) for NSD events 
at a/s = 0.9 and 2.36 TeV. The increase of (28.4 ± 1.4 ± 2.6)% from 0.9 to 2.36 TeV is signifi- 
cantly larger than the 18.5% (14.5%) increase predicted by the pythia (phojet) model tunes 
used in this analysis. The collision energy dependence of the measured dN^h/ drj\^^Q is shown 
in Fig. 03, which includes data from the NAL Bubble Chamber [25J, the ISR [26J, and UAl [22J, 
UA5 [24J, CDF EZI, STAR ||28l, PHOBOS ||29l and ALICE [23J. The dN^h/d}] measurement 
reported here is consistent with the previously observed trend. 



7 Summary 

Inclusive measurements of charged-hadron densities, cfN^h/ dpj and dN^h/ df], have been pre- 
sented based on the first pp collisions recorded at = 0.9 and 2.36 TeV by the CMS ex- 
periment during LHC commissioning in December 2009. The numerical values of the data 
presented in this paper can be found in Ref. |30|. For NSD interactions, the average charged- 
hadron transverse momentum has been measured to be 0.46 ± 0.01 (stat.) ± 0.01 (syst.) GeV/ c at 
0.9 TeV and 0.50 ± 0.01 (stat.) ± 0.01 (syst.) GeV/c at 2.36 TeV. The three reconstruction methods 



15 



employed for the dN^h/ dt] measurement have yielded consistent results, demonstrating the ex- 
cellent performance and detailed understanding of the CMS tracker. The pseudorapidity den- 
sity in the central region, dN^h/ d?/||^|<o.5/ has been measured to be 3.48 ±0.02 (stat.) ± 0.13 (syst.) 
at 0.9 TeV and 4.47 ± 0.04 (stat.) ± 0.16 (syst.) at 2.36 TeV. The results at 0.9 TeV have been 
found to be in agreement with previous measurements in pp and pp collisions. With the new 
measurements at 2.36 TeV, which show a steeper-than-expected increase of charged-hadron 
multiplicity density with collision energy, the study of particle production in pp collisions has 
been extended into a new energy regime. 




Figure 7: (a) Energy dependence of the average transverse momentum of charged hadrons. 
The CMS data points are evaluated for the range \f]\ < 2.4. Data of other experiments are 
taken from Refs. I|22l |3T1 - I33l . The curve shows the fit to the data points of the form (px) = 
0.425 - 0.0197 ln(s) + 0.00156 In^(s) with (pj) in GeV/c and s in GeV^. The error bars on the 
CMS data points include systematic uncertainties, (b) Charged-hadron pseudorapidity density 
in the central region as a function of centre-of-mass energy in pp and pp collisions including 
lower energy data from Refs. ||22i - |29l , together with various empirical parameterizations fit 
to the data corresponding to the inelastic (solid and dotted curves with open symbols) and to 
the NSD (dashed curve with solid symbols) event selection. The error bars indicate systematic 
uncertainties, when available. 
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